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NATIOKAL ADVISORY COMMITTEE FOR AERONAUTICS 
RESEARCH MEMORANDUM 

APPLICATION OF WING -BODY THEORY TO DRAG REDUCTION 
AT LOW SUPERSONIC SPEEDS 

By Barrett S. Baldwin, Jr., and Roller t R. Dickey 

SUMMARY 


A method is developed for extending to hl^er Mach numbers the region 
of low drag attainable for wing-body combinations by the use of the tran- 
sonic area rule. It is found that to a good approximation, the drag 
depends only upon the longitudinal distributions of area and moinents of 
area about the vertical plane of symmetry parallel to the free-stoeam 
direction. The essential requirement of the method la that the longi- 
tudinal developnent of the moments of area be smooth and gradual. 

Results of an experimental investigation conducted in the Ames 2- by 
2-foot transonic wind tunnel to test the theory are presented. The 
results in essence confirm the predictions of the theory in that the zero- 
lift wave drag of a wing-body configuration over a range of low supersonic 
Mach numbers as well as at sonic speed is reduced when auxiliary bodies 
are mounted on the wing. 


INTRODUCTION 


R. T. Jones has expressed the theory of wing-body wave drag at super- 
sonic speeds in a form which Illustrates the dependence of the drag upon 
the longitudinal distributions of the cross-sectional areas of the ccxnplete 
configuration intercepted by planes inclined at the Mach angle of the flow 
( see ref . l) . The derivation contains as a special case for a Mach number 
of one the transonic area rule Introduced by V/hitcomb (ref. 2) wherein the 
intercepting planes are normal to the longitudinal axis of the configura- 
tion. It was concluded In reference 1 that the modification of a wing-body 
combination in accordance with the transonic area rule would generally be 
expected to result in drag reductions at near sonic speeds; however, it 
was pointed out that, at hl^er supersonic Mach numhers, this modification 
would sometimes result in drags greater than that of the original config- 
uration. In reference 1 a method for contouring the fuselage of a wing- 
body combination was presented which achieved drag reductions at particular 



2 


NACA RM A54J19 



supersonic design Mach numbers, but only at the expense of Increasing the 
sonic-speed drag compared with that of the corresponding transonic-area- 
rule configuration. Thus, It appears that the methods for minimizing wave 
drag described in references 1 and 2, which are both based on the longi- 
tudinal distributions of cross-sectional area, are effective for only a 
limited Mach number range. 

In the present paper, Jones’ generalized zero-lift wave -drag formula 
is re-examined in an attempt to develop a method for minimizing the wave 
drag of a wing-body combination over a wider range of Mach numbers. 


ANALYSIS 

Calculation of Zero-Lift Wave Drag 


It has been pointed out in reference 1 that the transonic area rule 
was predicted by the linear theory, but was discounted because basic 
assumptions of the theory are violated in this application. It has been 
suggested by Jones that other predictions of the linear theory which may 
have been overlooked should be systematized and investigated experimen- 
tally. Thus, in the present analysis a possibly unwarranted emphasis is 
placed on the formal predictions of the linear theory at Mach numbers neetr 
one. 


In reference 3 methods are given for calculating the aerodynamic 
forces on airplane configurations utilizing very few assumptions other 
than those needed for linearization. An additional approximation is 
employed In reference 1 to relate the supersonic zero -lift wave drag of 
a wing-body combination to the drags of a series of equivalent bodies of 
revolution each of which is determined from the cross-sectional areas 
intercepted on the configuration by a set of parallel Mach planes. The 
result of reference 1 coincides with the more exact result of reference 3 
at sonic speed, and the deviation with increasing Mach number is expected 
to be small in a limited range of Manh numbers as long as the configura- 
tion is a conventional monoplane type. 

In the interest of obtaining a result in terms of familiEir geometric 
concepts and to facilitate calculations, the method of reference 1, tenned 
the "Mach plane method," will be employed here. This approximate theory 
greatly simplifies the discussion of gross effects of variations in the 
design of wing-body combinations. 

As a preface to the development of the method of the present paper 
for calculating drag, the Mach plane method will be briefly reviewed. 
Symbols are defined in Appendix A. 
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Mach plane method .- Consider a wing-hody ccmbinatlon such as shewn 
In sketch ^.a) • Let x he the coordinate In the free-stream direction, 
y the spanwise coordi- 
nate, and 2 the remain- z 

ing Cartesian coordinate 
in the thickness direc- 
tion, with the origin at 
the center of the hody. 

A Mach plane can he 
defined as a plane with 
its normal at an angle 
of tan-i(l/p) to the x 
axis. Let (x',p,cp) denote 
the Mach plane which inter- 
sects the X axis at x* 
and has the projection of 
its normal on the yz 
plane at an angle q) to 
the y axis . Let 
s(x',P»cp) 136 the area of Sketch (a) 

the projection on the yz 

plane of the cross-sectional area intercepted on the configuration hy the 
Mach plane (x*,P,cp). Then the drag of the configuration is the average 
with respect to cp of the drags of the equivalent hodies of revolution 
defined hy the area distrlhutions S(x*,p,cp). 

A method Introduced in reference It is used in reference 1 to evaluate 
the drag of each equivalent hody of revolution. The variable 0 is 
defined hy the relation 



x' 


— cos 0 
2 


CD 


where I is the length of the equlval^at hody. Then a set of quantities 
An(p,cp) are defined as the coefficients of sin n0 in a Fourier series 

expansion of Consequently the Aa(P,cp) can he determined 

from the relation 


2 

n 



sln(nS)d0 

Sx* 


An(P,<P) 


( 2 ) 
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Finally, the drag of the configuration is given by 




Within the framework of the linear theory this result is valid only for 
equivalent bodies of revolution witii no discontinuities in the gradients 
of the area distributions. 

It should be noted that unless all parts of the configuration lie 
between the nose Mach cone and the forward Mach cone from the tail, the 
equivalent body length, Z, will be greater than the actual body length 
in some cases. However, by consideration of streamwise body extensions 
of vanishingly small cross-sectional area, it can be seen that a constant 
value of I equal to or greater than the length of the longest equivalent 
body can be used in equation (l) . 

Series-expansion method .- In this section the Fourier series coef- 
flclents defined in equation ( 2 ) will each be expanded in a finite series 
so that the drag formula can be expressed as a power series in powers of 
3. This manipulation leads to an expression of the drag in terms of a 
convenient set of geometric parameters which were not apparent in the Mach 
plane method. 

By the use of equation ( l) , equation ( 2) can be written an 


An(3,(p) = - II f 


-Z/2 

or after a partial integration 


sin(ne) 
sin 6 


/' /“ Z / 2 _ 

“^nCe^cp) =(t) - / S(x',3,cp) 


sin nB 
sin 0 


provided that and S(x’,p,q>) are zero at the nose and tail.^ 

ox’ 


^For a practical configuration where the area distribution is not 
zero at the tall, the distribution of a Karman ogive having a base area 
and length equal to that of the configuration under consideration can be 
subtracted from S(x’,3,<p) so that the resulting equivalent -body area 
distributions will be zero at the nose and tail. The drag due to the 
part removed can then be calculated by means of equation (2) rather than 
equation ( 4) . The choice of a Karmm area distribution insures that there 
will be no interaction drag from the part removed as long as all configu- 
ration parts lie within the Mach cone from the body nose eind within the 
forward Mach cone from the tall, as can be seen by the use of equa- 
tions ( 2 ) and ( 3 ) - 


V 
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The second factor of the integrand of equation {k) 

X* 

in a finite series of powers of given by- 


can be expanded 


n-2 




sin ne 
sin 6 




n> 2 


m=o 


■where 


Jnm - 


Me; 


for even values 
of n - m 


zero for odd values of n - m 


( 5 ) 


C6) 


Substituting equation (5) into ( 1 ^) and interchanging the order of 
summation and integration yields 


An(^,9) = ^ I J S(x«,p,cp)(^^j dx« , n>2 (t) 

m=«D -l/s 

At this point, a more explicit expression for S(x’,p,cp) in terms of 
the configuration geometry is needed for substitution in equation (t) • 

Let t(x,y) be the thickness distribution of the configuration including 
that of the body or bodies. It -will be assvuaed that the distance of all 
parts from the xy plane is small enough that S(x’,p,cp) can be approxi- 
mated by 


ryaCx' ,P,cp) 

S(x',p,cp) = J t(x' + 3y cos cp,y)dy (8) 

yi(x',p,(p) 

■where yi(x*,p,q)) and yaCx'^P#^) define the edges of the configuration 
intercepted by the Mach plane (x’,3,9). Equation (8) represents a planar 
approximation. Analogous expressions not involving this ap-proximation 
can be found and exploited, but only the planar case -will be discussed in 
S this report. 







Subatitutioa of equation (8) Into equation ( 7 ) yields 


^ / \2 r^/2 oy£(x’,p,9) , 

An(p><P) = 2^ S. J J t(x' + py cos Cp,y)dy (9) 

m=0 x'--7/a y=yi(x^p,(p) 


If it is understood that t(x,y) is zero at points off the configuration rather than an analytic 
continuation of its form at points on the configuration, the integrations with respect to x’ and 
y can both be taken from -<» to a> and the order of integration interchanged. In addition, with 
the Eubstitution of x = x’ + py cos <p, equation (9) beccmes 


An{p,ip) = ^ (ej I J t(x,y) (\ 


X - py COB 


y=-oo x=' 


djc dy, n^a 


The quantity in parentheses which is raised to the power m can be expanded by the blncraial 
I theorem into 




where 


^ m' 


Then substituting equation (ll) into equation (lO) yields 


MM -I if f I (-1)' (^J"' ay 


y-~oo x=*-» p=o 


t -K 
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Upon interchange of orders of integration and snunnation, shortening of the interval of integration, 
and arbitrary grouping, this becomes 


A«(P,<P) = f ‘T* (-1)' “p (fj I (i^)“ 

m-o pBO , x=-2/e Ly»-ya(x) J 

where y„(x) and y 4 (x) define the two edges of the configuration. 


dx (13) 


The quantity in brackets in equation {13) can be identified as the longitudinal distribution 
of the pth moment of area of the configuration. This Indicates that the drag of the configura- 
tion can be expressed entirely in terras of moment distributions (including the area dlstrlbutlan 
which correspands to p ^ o) . 

The moment distributions can be defined as 

ry^Cx) 

Mp(x) « J t(x,y)yPdy (l 4 ) 

y=-y0(x) 


Substituting this In equation (13) yields 


n-a m 


An(p,q)) = ^ bum ^ (-1)^ cj | J Mp(x)^Y^^ ^ pPcosPip 


dso 

or interchanging the order of summatlonB yields 


An(P,«p) - ^ (-l)^PcosP(p ban I (f)^^ ^ 

pp=o Lm=p -l/s _ 


as the desired expansion of the An(p,q>)'s in powers of p. 
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It can be seen in equation (l4) that if the configuration has span- 
wise symmetry, the odd moment dlstrlhutlons will he Identically zero, and 
the terms of equation (l6) resulting from odd values of p will he zero. 
With streamwlse symmetry of the moment dlstrlhutlons In addition, the odd 
values of n and m would he eliminated. 

In the process of substituting equation (l6) into Jones' drag equa- 
tion, it Is convenient to define several new symbols. 

Let ■ ■■■ 


n -2 


m=p 


In' o / \m-p 

Inp = X C- I (f) f Mp(x)(^) dx ( 17 ) 

-2/e 


so that 


n -2 


An(3,cp) = ^ (-l)’^p^cos^cp Lnp 


p=o 

Then can be written as 


(18) 


n -2 n -2 


[A: 




^o(Pi+Pe) (q,)|3(Pi+Pa) 


.cos 


Pi=o Pa=o 


or 


f [An(3,(p)]^ = 


2n-4 

I 

q=o 

q even 


q -2 

2 


+2 ^ Lnpi Ln(q-pi) 


Pi=q-n+2 


cos‘3<p3‘1 


where the odd values of q are omitted because the terms resulting from 
such values would not contribute to the drag. 

In addition to Lnp of equation (l?)# there are several other 
quantities depending upon Lnp which are notatlonal aide. Let 



( 20 ) 



NACA RM A 5 lWri 9 


9 



and 


Then 


Slz£ 


s 



Pi=q^-n+a 


( 21 ) 


an.-'t 

[An(P,<p)l^ = 

q=o 

q even 



+ I^^)cosl<pp<l 


Substituting this In equation (3) yields 


( 22 ) 


an — 4 


D 


2 lt , 

= ^ f Z ^^aq + Inq)cos<lcppiajtp 


n =2 


q=o 

q even 


or, if it is assumed that the series may be integrated term-by-term, this 
becomes 


D = 


cPV^ 

8 


00 2n-4 


^ n ^ (Nnci + Inq) JqP ' 


(23) 


n=2 q=o 

q even 


where 


^q. 


= i f 

2jt J 


2 It 


coslcp dtp = 


q.£ 


(t>r 


, q even 


(24) 


Interchanging the order of summations in equation (23) results in 


D = 


itpV^ 

8 


00 

Z Jq.P‘ 


q=o 

q even 


00 

Y n(Nnq 


+ Inq) 


n= 3 il 

a 


(25) 


DEN^IAL 
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Let IJq te the quantity derived from Nnq lay 

00 

Nq = ^ n Nr 


’nq 


n= 




and define Iq as 


= z 


n X 


■nq 


n=3^ 

2 

Then substituting these in equation ( 25 ) yields 

00 

D = ^ Jq(N(^ + Iq)^'l 


<J^5sO 

q, even 


( 26 ) 


( 27 ) 


( 28 ) 


as the desired expansion of the supersonic zero-lift drag formula in 
powers of 3 . 

From the foregoing, it is seen that each Nq depends only upon the 
longitudinal distribution of the moment of area of order q/2, whereas 
each Iq depends on all moment distributions of order zero to [(q/s) - 1]. 
Thus, each Nq represents a contribution to the drag from the moment dis- 
tribution of order c[/2 alone, and each Iq represents a contribution 
resulting from the interaction of the first [(q/2) - 1] moment distribu- 
tions . 

Although the question of convergence of the series of the foregoing 
analysis has not been investigated in detail, several observations and 
practical hints for calculation can be offered. 

The values of each An(3,'p) obtained by the series-expansion method 
are identical to those obtained by the Mach plane method. Therefore, if 
it is assumed that the drag of a configuration can be calculated with stif- 
ficlent accuracy by using the first N terms of the Mach plane method, it 
follows that the interchange in order of summations by which equation ( 25 ) 
is derived from equation ( 23 ) is valid for these terms. 

In the Fourier series analysis of 9s(x ' ,3 ,cp)/5x' in the Mach plane 
method, it is evident that the higher harmonics will be suppressed if the 
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smallest allowed value of I is used rather than a large value. Conse- 
quently, the convergence of hoth the Mach plane and series -expansion 
methods is hest when the smallest allowed value of 1 is used in equa- 
tions (l) and (it) • In either method the number of terms required to 
obtain the major part of the drag can be held to a minimum, and the 
mathematical calculations thereby facilitated, by dividing the configura- 
tion under consideration into a short part and a long part. The quanti- 
ties .^CPjCp) of the complete configuration are the sums of the corre- 
sponding quantities of the separate short and long parts as given by the 
relation 

An(Pjtp) = ASn(Pj'P) + AltlCP^t) (29) 

where the subscript S is used to denote short part and L long part. 
Then [An(p,qp)]^ is given by 

[An(3,cp)l^ = [Asn(P,cp)l^ + 2Asn(P,cp)ALnO,cp) + [ALn(3,cp)]^ (30) 

When this expression is substituted in equation ( 3 ) it is seen that the 
[Asn(Pj<p)]^ will yield the drag of the short part alone, [ALn(P,«p)]^ 
the drag of the long part alone, and 2Asn(^;tp)ALn(^,q)) the interaction 
between the short and long parts. Then a smaller value of 1 can be 
used to calculate the drag of the short part alone. Although no reduc- 
tion in the value of I is possible for the other two parts of the drag, 
the convergence is Improved because of the absence of the high harmonic 
content of the short part from AT.n(P.tp) • 

The number of terms to be Included in the drag formula will depend 
upon the relative importance of accuracy and simplicity. In the search 
for low-drag design criteria, a very small number of terms might be 
appropriate. As an example of the meaning of this remark, it can be seen 
that only two terms need be considered to arrive at the supposition that 
the Sears -Haack area distribution is an optimum for given length and 
volume . Only one teirm is needed to conclude that the fineness ratio of 
a body should be as large as possible when the pressure drag alone is 
considered. 


A Method for Reducing Drag 


The general problem which will be considered in this section is that 
of designing a wing -body combination with low drag in a range of super- 
sonic Mach numbers when certain basic parameters, such as total volume, 
are specified. In order to obtain definite answers, many additional 
parameters such as those involved in the specification of the wing plan 
form must be assigned arbitrarily. For example, the wing can be chosen 
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arbitrarily, and the optimum body shape for minimum zero-lift wave drag 
at a specified Mach number can be fotmd by methods described In refer- 
ence 1. 

If the number of parameters affecting the drag Is small, as Is the 
case when the transonic area i*ule Is valid, the general problem of drsig 
minimization la greatly simplified. The series-expansion drag formula of 
the present analysis is also expressed in terms of a small number of 
parameters if the higher powers of g can be neglected or if the stmima- 
tlons over n can be cut off at a small number. Because of the resulting 
simplification it has been foiind that the minimization procedure employed 
in reference 4 can be used to design an optimum wing-body combination 
with minimum drag at a Mach number of one and minimum drag rise at low 
supersonic speeds. As a first application of ‘ the foregoing analysis this 
procedure will be, described and exploited. The validity of the result Is 
subject to some .question because the basic assumptions of the linear theory 
are violated at Mach numbers near one. However, the results are of inter- 
est in the absence of a method for applying a more exact theory. 

For the ordinary case of spanwiee symmetry, equation (28) can be 
written as 


(31) 


In the speed and aspect-ratio range where p® and higher powers of 
P can be neglected, the drag depends only upon the area distribution, 
since this is the only feature of the geometry affecting Nq. It can be 
assumed that the geometry will be such that the hl^er powers of p 
should be taken into account successively as the value of p is increased. 
Then in the speed range where powers of P greater than two can be neg- 
lected, the drag depends only on the area distribution and the second- 
moment -of -area distribution, since these two determine the value of la. 

As successively higher powers of p are taken into account, correspond- 
ingly higher ordered moment distributions are Involved. 

The following procedure is proposed for reducing the wave drag of a 
wing-body configuration: 

1. Minimize the drag at a Mach number of one by exclusive attention 
to the area distribution. 

. Minimize the drag at slightly hi^er Mach numbers by exclusive 
attention to the second-moment-of-area distribution excluding 
any changes which would alter the area distribution. 


D = Jiez! 
8 


Wo + I + 
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3 . Continue to find the optimum hi^er ordered moment distributions 
successively without disturbing the lower distributions in 
their previously derived optimum forms. 

This procedtire leads to a unique set of moment distributions. 

Substituting equation ( 20 ) into (a6) with q_ = O yields 


W 

= ^ (bno) 


Frcan equation (l^) with p = 0, the volume of the configuration, Vq, is 
identified as 

pi /z 

Vo = J Mo(x)dx (33) 

- 1/2 

By the use of equation (17) with n = 2 and p = 0, t^e first term of 
equation ( 32 ) is found to be proportional to (Vq/Z^ . Then, since a 1 1 
the terms of equation (32) are positive, Nq is a mlnlinum for given 
values of Vo and I if 


Luo = 0 for n / 2 (34) 

To satisfy equation (34) and similar equations which will occur, a 
rearrangement of equation (I 7 ) is needed. When the integral and sum of 
this equation are reversed to obtain 



the quantity in brackets can be identified as 
with the aid of eqtiations ( 5 ) and (l2) so that 



sin n9 
sin 0 
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Then after (p+l) partial Integrations this becomee 


2 2 Mp(x) 

^P-?T / ,_ 


-Z/2 


sin n 6 
sin 6 


dx 


(35) 



for all integer values of q frcan zero to p. Equation ( 36 ) is satis- 
fied if Mp(x) and the first p derivatives of Mp(x) are zero at 
X = ±1/2. 


It follows from equation (35) with p = 0 that equation (34) is 
satisfied by 


^ - Wip(ae) - Ji - (j^7 

By Integration and use of the fact that the conflginration does not extend 
beyond x = ±l/2, this becomes 

Mo(x) = i I L 20 [1 - 

For purposes of evaluating ttie drag, this area distribution can be put in 
the form 




Mo(x) = Mo(o) 


(39) 
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trtiere Mo(o) is the maximum value of the dlstrlhution . Then 



m 


<and 


D 



( 1 ^ 1 ) 


Is obtained as the drag of the optimom configuration in the speed range 
■where 0 can he neglected. Equations (39) ani (^l) are In agreement 
■with the results of references Ij- and 5* 

By suhstitution of eguntion ( 3lt-) Into equations ( 2l) and ( 27") , it 
is found that Ig and I 4 ■will he zero for a configuration with the opti- 
mum area distribution and all other Iq's ■will he Independent of the 
area distrihution. In that case equation ( 3 I) becomes 




~Mq(o)1 2 

I 


+ 




+ 


^(Ne + Is)P® + 0(pi°)} 


(42) 


In minimizing it can he assumed "that the second -moment distri- 

hution of the hody is negligible so that I can he replaced by Ig, the 
■wing length. More exactly, Ig Is the length of the projection of the 
■wing on the hody axis. Substituting equation (20) into ( 26 ) Vith q = 4, 
yields 

00 

n=4 



The second-manent volume of the configuration can he defined as 


Vg = 



-W2 


Mg( x) dx 


(44) 
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By the use of equation (ij) with p = 2 , the first term of equation ( 43 ) 
is found to he proportional to Then since all terms of equa- 

tion (43) are positive, N4 is a minlnum for given values of Vs and 2 a 
if 

Lq 2 = 0 for n ^ 4 (45) 

It follows from equation (35) with p = 2 that this requirement is met hy 
setting 


d^Ms( x) 
dx^ 


= L42sin(45) 


(^) 


By triple integration and use of the sufficient requirement of equa- 
tion (36) that Ma(x) and the first two derivatives of MaCx) he zero at 
X = ±22/2, this becomes . _ . . ...1_ 


Ma(x) 



( 4 t) 


or 


M^(x) = 


(M) 


is obtained as the optimum second-moment distribution. Then 

MaCo) 


L42 = 105 


(^2)® 


( 49 ) 


and by substitution of this in equations ( 42 ) and ( 43 ), the expression 


D = 


XPV2 


{18 


Mo(o) 


+ 16537.5 


Mg(o)~ |gp4 

UZa)®J 


(50) 


is obtained as the drag of a wing-body combination for which the distri- 
butions of area and second mcment of area are optimum when powers of p 
greater than four are neglected. 

The foregoing process can be continued indefinitely until a complete 
set of optimum moment distributions is obtained. The restilts are 
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l^(x) = Mp(o) 


b - (^)T 


+(s/e) 


(51) 


where 


lo = I 

Ip = la for p >2 '-7 


arfl the drag of such a configuration is given hy 


D = 


itpyg 

8 


00 

I 


p=o 

p even 


[ ( gp) '. ] H gp+l) ^ Mp(o) 

(p!)7b+2) 


p2P 


( 52 ) 


Since all the terms of equation (52) are positive, this drag a.nrt the 
corresponding drag coefficient must increase monotonically with increasing 
Mach numher. Also since all the terms are positive it can he seen from 
the first two terms that the drag will he very large at 




Mo(o)_^ 

Ma(o)l 


(53) 


Equation (53) caJi he used to estimate the upper limit of the Mach numher 
range of applicahillty of the foregoing low supersonic technique for drag 
reduction . 


It is interesting to note that a configuration designed according to 
equation ( 5 I) would have large drag at the hlgsher Mach numhers as a con- 
sequence of eliminating the drag components due to interactions of the 
moment distrihutlons . Although the drag due to each moment distribution 
alone must he positive, the interaction drags can he negative. The inter- 
action drags which would he beneficial in a given Mach numher increment 
were eliminated in the process of minimizing the drag at lower Mach num- 
hers . 


FlUENTlAL 
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APPLICATIOM'S ABD DISCUSSION 
Design of Configuration for Low Drag 


In the "analysis" section, it has heen shown that the wave -drag 
formula can he expanded in a power series of the form 


D = Sq + + S4S4 + . . . (54) 

where the constants so forth, are Independent of Mach number 

and are determined only from the geometiry of the configuration. In this 
section, the physical significance of these constants will he discussed and 
some examples presented of the practical means available for minimizing the 
constants in order to reduce the drag of wing-body combinations at low 
supersonic speeds as well as at Mach number one. 

The details of the procediire for deriving the constants frcan the 
geometry using a planar approximation are contained in the "ANALYSIS" 
section and in Appendix B. Equations (ll) , (.SO), (21), (26), (27), and 
(28) indicate that, in general, ao depends only on the sirea distribution, 
ag depends upon the second-moment- of-area distribution as well as on the 
area distribution, and a 4 depends on the fourth-moment- of- area distribu- 
tion in addition to the previous two distributions. Hence it is seen that 
the transonic area rule is valid in the speed range where p is small so 
that all terms except the first in equation (54) can be neglected. Fur- 
thermore, it is expected that an the Mach number is increased, starting 
from one, all except the first few terms should remain negligible in a 
range of low supersonic Mach n\imber 3 so that at these speeds the drag 
should depend only upon the area distribution and the second-mranent- of- 
area distribution. 

If the configuration area distribution is made an optimum for minimum 
drag at a Mach number of one by the use of the transonic area rule, the 
determination of and 34 is simplified. In that case ag is zero and 
a 4 depends only on the second-moment-of-area distribution. Consequently, 
if the second-moment-of-area distribution can be varied without changing 
the area distribution, the drag at low supersonic speeds can be minimi zed 
with respect to . such, variations without increasing the sonic speed drag 
which depends only on the area distribution. In order to see that the 
second-moment distribution actually can be varied without changing the 
area distribution, definitions of these distributions are needed. 

The area distribution is given by 
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■wtiere t(x,y) is the thickness dis trihut ion of the configuration incliod- 
ing wing and hody, and Y(x) is the value of y at the edge of the plan 
form. In this definition spanwise symmetry of the configuration has "been 
assumed. 


The second-moment distribution can he approximated hy 


MsCx) = 2 



y2t(x,y)dy 


(56) 


where R(x) is the value of y at the wing-hody juncture. Here the 
second-moment dlstrihution of the hody is neglected hecatise of the small 
values of y at the hody compared with those on the wing. Since the hody 
moments are negligihle, changes in the hody shape will vary the configura- 
tion area dlstrihution without altering the second-moment distribution. 
Conversely, the second-moment dlstrihution can he altered while holding 
the area dlstrihution fixed hy varying the wing geometry and the hody 
shape at the same time. 

If the area and second-moment distrihutions are made optimum, the 
drag can still he varied hy altering the magnitudes of these distrihutions, 
as can he seen in the drag formula for such a configuration given hy 


D = 9 jtpv2 

4 


^(°) + 3305 itpV^ ^g(°) + o(p®) 

_ ^ J 16 Ll T-s) 


(57) 


where Mq(o) Is the maximum value of the area dlstrihution, 1 is the hody 
length, Ma(o) is the maximtan value of the second-moment distribution, and 
Z 2 is the wing length. 

Neglecting powers of p greater than foirr in equation (57) provides 
an insist Into the requirements for reducing the pressure drag at low 
supersonic speeds. The area-rule requirement that the ratio Mo(o)/l he 
small indicates that the fineness ratio of the hody should he large as 
previously noted. Examination of the quantity M 2 (o)/Z 2 ® leads to the 
conclusion that not only should the thickness ratio of the wing he small, 
hut also the thickness ratio should taper to a minimum at the wing tips, 
and the ratio of effective strearawise length to span of the wing should 
he large. At hl^er Mach ntunhers, where the higher powers of p cannot 
he neglected, these conclusions would not apply. 

The ratio of effective length to span of a wing can he Increased in 
several different ways, for example, hy extending the wing chord. However, 
large frictional drag penalties are usually associated with the Increased 
surface area accompanying such changes. Another possible method of 
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Increasing the effective-length-to-span ratio of the wing is hy the siddl- 
tion of auxiliary bodies of revolution moimted on the wing. This method 
has the advantage of relatively small increased surface area and attendant 
friction dra^. The bodies of revolution are particularly attractive in 
the case of wing-body configurations of relatively large wing span where 
application of the transonic area rule could be expected to produce drag 
penalties in the low supersonic, speed region. 

As an illustration of the application of the drag-reduction proce- 
diire, hereinafter referred to as the ”moment-of-area rule" as distin- 
guished from Whitcomb's area rule, consider the wing-body combination 
shown in the upper part of sketch ("h) • This configuration consists of a 

actual distribution 

-------- optimum distribution 


w. 



combination Sears -Haack-Karm^ ogive body of fineness ratio 11, and an 
elliptic -plan -form wing of aspect ratio 2.0 with circular-arc sections 
and 5-percent maximum thickness ratio. The distributions of area and 
second moment of area for this basic configuration are also shown in the 
sketch. The shapes of these distribution curves are not conducive to low 
drag in that the area distribution has a bump at the location of the wing 
and the moment -of -area distribution is short and has steep slopes. With 
the total volijune fixed, the optimum shapes of the distribution curves (as 
defined by eq,. (51)) Bxe shown by the dashed lines. ^ The desired dlstribu- 
tion of the second moment of area can be obtained by utilizing auxiliary 
^The optimvim second-moment distribution is not a function of Mach 
number because it is derived essentially by minimizing the derivative 
dCp/dP'* evaluated at a Mach number of one. 
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(NITDENTIAL 


todies , of revolution mounted on the wing as shown in the second part of 
sketch (h) > The arhitrarily chosen spanwise location of the auxiliary 
todies determines their size in that small todies at an outtoard position 
can produce the same second-moment dlstritution as larger todies at .an 
Intoard position. It is evident that in order to prevent an increase in 
the maximum value of the second moment of area the auxiliary todies must 
te waisted in the vicinity of the maximum thickness of the wing. The 
area dlstritution may te made optimum ty reshaping the tody to satisfy 
the requirements of the transonic area rule after the aixxlllary todies 
have teen added. 

In discussing the effects of modifications it is convenient to iso- 
late portions of the drag which will not te affected ty the modifications 
\mder consideration. Considering presstire drag only, the quantity of pri- 
maiy interest is the additional pressure drag caused ty all additions to 
and alterations of the original tody alone. The wing and auxiliary todies 
are considered to te additions while the reshaping of the tody is an 
alteration. Another reason for isolating this additional pressure drag 
(ACp) is that the tasic assumptions of the linear theory used to calculate 
ACp for configurations with the transonic-area-rule modification may not 
te violated, although the aasxmiptions are violated at Mach rrumters near 
one for the tody alone (see ref. 3 ). 

The additional pressure drag as Just defined is ottained ty calcu- 
lating the drag of a configuration consisting of the wing, the auxiliary 
todies, and the tody cutout. The tody cutout is taken to te a negative 
area dlstritution located at the position of the tody surface. The calcu- 
lated values of ACp for the unmodified configuration, and the configura- 
tion modified according to the mcsnent-of -area inile are shown in sketch ( c) . 



Sketch (c) 



22 


NACA RM A54J19 



For comparison, the calculated values of ACp for a configuration modi- 
fied according to the transonic area rule only are also shown. A descrip- 
tion of the methods employed in the calculations appears in Appendixes B 
and C. The drag coefficient is based on the total ving area including the 
peirt of the wing hidden inside the body. It is apparent that the addition 
of the auxiliary bodies to mahe the second-moment-of-area distribution an 
optimum results in large theoretical drag reductions at low supersonic 
speeds . It is to be expected that the actual dretg reductions will be 
somewhat less than those predicted because of the effect of friction drag 
not taken into account by the theory. , 

Experiment 

In order to obtain an experimental check of the theoretical predic- 
tions, models of the configurations under consideration were constructed 
.and tested in the Ames 2- by 2-foot transonic wind tunnel at a Reynolds 
niimber of I .9 million based on the wing root chord. 

The experimentally measured values of the total drag coefficient at 
zero lift for the configurations are shown in sketch (d) . 

O unmodified wing-body configuration 

□ area-rule modification 

O moment-of-area-rule modification 


^ unmodified body alone 



M 

Sketch (d) 
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These results canfirm the predictions of the theory In that the drag 
of the configuration modified "by the -use of atixillary hodies is nearly as 
iow as that for the configuration modified hy the area rule alone at a 
Mach numher of 1.0 and is less at higher Mach numhers. The drag of the 
configuration ■with auxiliary "bodies is greater at su"bsonic speeds than 
that of the other configurations "because of the Larger surface area. 

The drag of the unmodified hody alone is also shown in sketch (d) 
together with the predicted supersonic value which is plotted as an incre- 
ment above the experimental value of drag at low speeds. The poor agree- 
ment is considered to result from a reduction of skin friction at super- 
sonic speed due to an Increase in the extent of the laminar flow since at 
supersonic speeds the extended regions of falling pressure are conducive 
to delay of transition to tur"faulent flow. 


The Incremental drag rises with the body drag excluded from the 
experimental and theoretical values show much better agreement as can be 
seen in sketch (e) . 





area -rule 
modification 


unmodified 


moment -of -area- 
rule modification 


Sketch (e) 

The ZiCp values were estimated from the experimental results by sub- 
tracting the subsonic drag of the ccaaf iguratlon as well as the drag rise 
of the original body alone from the total drag of the configuration. This 
operation can be expressed as: 

= ^D-fcotal “ (^^®to-^al^g^ M=o*s ” t*^"’^tiody alone 



dy alone 




at M=o. 


6 




? nv^!lfr rT5T:^ 
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or 


(^*^^total ~ ^^ody alone^ “ ^*^®total “ *"^ody alone^^^ ^ 

For the auxiliary bodies employed in the foregoing experimental 
investigation, the assiimption that powers of p greater than fotir can 
he neglected in the drag formula is violated at Mach numbers above 1.1. 
Hence it is expected that modifications more effective in drag reduction 
at the hi^er Mach numbers can be found if this assumption is not used, 
■fhe series -expansion method for evaluating the drag can be used to design 
auxiliary bodies which will minimize the pressure drag at a specified 
sttpersonic Mach number if the higher values of n are neglected rather 
than the higher powers of p. The result would correspond to the fuse- 
lage modification for minimum drag at a specified supersonic Mach number 
described in reference 1. 


COWCLUDING REMARKS 


A basic method for estimating the first-order deviations of the drag 
of wing-body combinations from the values predicted by the transonic area 
rule has been derived. In a planar approximation it has been found that 
at Mach numbers above one the zero-lift wave drag depends on the distribu- 
tions of moments of area of the configuration about the vertical plane of 
symmetry parallel to the free-stream direction as well as on the area 
distribution. Thus the area rule can be supplemented by what ml^t be 
termed a moment -of -area rule for extending to higher Mach numbers the 
drag reductions associated with the use of the area rule at a Mach number 
of one . 

Just as Is the case for the area rule where the longitudinal develop- 
ment of area must be smooth and gradual to minimize the drag, so also. In 
application of the mcmient of area rule, the longitudinal development of 
the moments of area must be smooth and gradual. It has been found that 
at low supersonic speeds the moment-of-area distributions of order higher 
than the second are of secondary importance. Significant drag reductions 
can be obtained at these speeds by mounting bodies of revolution on the 
wing for the purpose of Improving the second -moment-of-area distribution. 
This point hae been verified by an experiment performed in the Ames 2- by 
2-foot transonic wind tunnel. 

An alternative way to visualize the mechanism of drag reduction by 
this means is to regard the auxiliary wing -mounted bodies as local 
pressure-field cancellation devices in the same sense that Jones and 
Whitcomb employ the contoured principal body or fuselage to counteract 
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the pressure field for the entire wing. From this point of view, it 
should he expected that the auxiliairy-hody modification would he most 
applicable to configurations embodying wings of relatively large span 
where the area-rule effects would be limited because of the large dis- 
tances of some of the wing elements frcm the fuselage. 

The concept of introducing auxiliary bodies along the wing span to 
effect decreases in wave drag promises to find important application for 
aircraft intended to carry external stores. For such aircraft, the pos- 
sibility exists of shaping the stores according to the moment of area 
rule so as to obtain drag reductions at transonic speeds with no friction 
penalty at lower speeds. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Oct. I 9 , 195^ 
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APPEHDIX A 


TABLE OF SYMBOLS 


8 -n 


An(P»<p) 


"brinj 


coefficients of in a power-series expansion of the 

drag 

coefficient of sin(n0) in a Fourier series expansion of 
9 s(x' ,$,cp) /^x' (See eq. ( 2 ).) 


coefficient of 



in a power-series expansion of 


5 sin(ngVsln(9) 



(See eqs. (5) and (6) .) 


h 

Cf 


D 


I 

N 


q^^inq’l 

qj^nqJ 



I 

I2 


M 

Mq(x) 

Mg(x) 


wing span 

coefficient in a hlnomlal expansion (See eqs. (ll) and’ 

( 12 ).) 

zero -lift wave drag 

quantities Involved in evaluation of the drag 
(See eqs. (20), (2I), (26), ( 2 T), and (28).) 

(See eqs. ( 2 ^ 4 -) and (28),) 

coefficient of ( S cos cp)^ in power-series expansion of 
Aji(p,<p) (See eqs. ( 17 ) and (I8).) 

hody length 

length of the longitudinal dlstrihutlon of second moment 
of area (l.e., the length of the projection of the wing 
on the X axis) 

free-stream Mach number 

longitudinal distribution of area of the configuration 

longitudinal distribution of mranent of inertia about the 
xz plane (also called the second -moment-of -area distri- 
bution) 
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Mp(o) 

Mp(x) 

R(x) 

S(x',p,cp) 

t(x,y) 

V 

Vp 

X 

y 

z 

(x',p,cp) 

yi(x’,p,cp)\ 

y2(x',p,cp)j' 

yaCx) jy^Cx) 
Y(x) 

p 

0 

p 



I 

I 


maximum veQ-ue of moment of area of order p 

longitudinal distribution of moment of area of order p 
(See eq^. (lit) .) 

dummy integers of siimmatlon 

Tralue of y at the wlng-hody juncture 

area of the projection on the yz plane of the cross- 
sectional area Intercepted on the configuration hy the 
Mach plane (x’,p,cp) 

thickness distribution of the configuration Including 
wings and "bodlee 

free-stream velocity 

pth moment volume (See eq^s. (33) and (^ 4 ).) 

Cartesian coordinate in the free-stream direction 

Cartesian coordinate in the spanwlse direction 

Cartesian coordinate in the thickness direction 

Mach plane which intersects the x axis at x’ and which 
has the projection of its normal on the yz plane at an 
angle tp to the y axis 

values of y at the points of intersection of the configu- 
• ration edges with the Mach plane (x*,p,(p) 

values of y at the edges of the configuration 

value of . y at the edge of a configuration with spanwlse 
symmetry 

speed parameter equal to 's/m® - 1 
variable related to x* by x* = ■g cos 0 

air density 

angle between the y axis and the projection on the yz 
plane of a normal to the Mach plane (x’,p,cp) 


«P 
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APPENDIX B 

DBiENSIONLESS EEAG FORMULA 


In order to take advantage of the decrease In the number of parame- 
ters resulting from similarity considerations, and to facilitate calcu- 
lations, the quantities defined in the "ANALYSIS" section can be made 
dimensionless. 

The moment -of -area distributions defined in equation (l 4 ) can be 
divided by a thickness to and the half-span b/2 raised to the (p + l)th 
power to obtain dimensionless moment distributions defined as 


Mp(x) = - 



tn(| 


y=-y 3 (x) 


^-Mp(x) /" (BI) 

Ha) 


similarly, double moments occiorrlng in equation (15) can be replaced by 

(B2) 


^ Mp(x)xi^dx 


(I) 




where c is the length of the projection of the wing on the x axis. 

With these definitions, a dimensionless version of the quantities 
Lnp of equation (l?) can be written as 


n-p-2 

Lnp = ^ Snpk » 


n - p - 2>0 


(B3) 


with the constants 


k=o 


Snpk 


defined as 


Snpk - Sbkp - 


(- 1 ) 


n-p-k-2 ( n+P+k' ^, _ 

2 V 2 y* 2^ 2P for even values 

j ^n-p-k-2^ , P- (n-p-k) 


zero otherwise 


(bL) 


i'liyiiN'r 
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Dimensiotiless Fourier coefficients can te defined as 

n-2 


A ]2 - 


I £ap(f -=?) 


2tol)/' 


p=o 
p 


and An® would te written as 


(B5) 


2 U-* 


An® = ^ Hnj cos cp^ 


where 


J=o 
S even 


d 


p=o 


(b6) 


(BT) 


p even p c-a/caa. t, 


Dimensionless drag components can he defined as 

\ZTt 


~ pzn 

Bn ~ — ^®dip 

2 tc 


(b8) 


so that 


where 


2 n-* 

Btt = n ^ Hnj Jj 

J=o 

J even 


1 ^ 1 ' 

Ja = — / cosJqp dip = iJ! 

^ 2it Jo ^ ^ r/A\ 1 


., J even 


20 (4) 


Then the drag Is given hy 




n-z 


(B9) 


(BIO) 


(Bll) 
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where 


q, = 1 PV^ (B12) 

The following is a suggested computing procedure: 

1. Evaluate the dimensionless double moments Mp^ for the configu- 

ration using equation (B 2 ) . 

2. Choose a value of Z equal to the length of_^the longest equiva- 

lent body of revolution and evaltiate the Lnn’s of equa- 
tion (B 3 ) , 

3. Evaluate the Hnj’s in equation (E7) and the Bn's in equa- 

tion (B9). 


The dimensionless drag components Dg, D 3 , D 4 , and so forth, should be 
evaluated separately so that the convergence with respect to n can be 
watched and computing errors found more easily. 


To determine which Mp^'s are needed for each value of n the sums 
can be written out proceeding in the opposite direction from that of the 
computing procedure. Thus 


D2 = 2H20 

B3 = 3H30 - ‘ ■ 



= 4 

H 40 

L. 

T 

+ |H42 

©■ 


1 1 

^ 


D 5 

= 5 

H 50 

1 ~ 

+ H 52 



©■; 


Be 

= 6 

Heo 

1 ~ 

+ |H62 


+ |h94 


^ (¥) 



(?J] 






etc 

• 






There will be 

(n - 

■ 1 ) terms in 

Dn for 

even 

n and (n - 2) 


odd n. 


\ (BI3) 

+ 


terms for 
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Equation (B3) indicates tliat the quantities L^p with p > n - 2 
are zero so that most of the terms in equations (b 6) are zero. Conse- 


quently , 
®ao 

^30 

= (Lao)^ 
= (Lao)^ 






A# 

A# 

A# 

A* 

«w 



^40 

•* 

II 

S 42 = 

2 L 40 L^2» 

H 4.4 = 

( L+a) ^ 



= (L5o)^ 

H 52 = 

2 L 50 Lsa, 

^54 = 

(Lsa) ^ 

)>(bi4) 

Heo 

= (Leo)^ 

tea = 

A« 

2I>ao Le 2 , 

II 

41 

(Lea) ^ ■*■ 2 Ij6o Ls4 


Hee 

^o> 

^ A* 

= 2Ls 2 Le^, 
etc. 

Hsa 

= (Le*)^ 


f ■ 



Erom equations (B3) and (bU) 


Lao = Moo 

Lao = ^01 

L40 = " ^o + 1^02 f L *2 = latao 

L50 = - la^ox + 32M0S > L52 = 9^21 

Leo = 3M00 ^ Mo2 + 8OM04 

Lea = - it8 M20 (^* 

= 8CM*o 





Lto^ etc. 

There will he — P 
terms for odd n. 


J 


terms in Lno for even values of n s.i 7 <^ 5 — ..E ~ ^ 

2 


‘OKFIDEIITIAL 
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Since only a small number of terms are Involved in Dn for n < 6, 
these fin's can he conveniently written out In terms of Epic’s as fol- 
lows: 



Because of the„small differences between large numbers Involved in 
the evaluation of ^Lnp for large values of n, it is probably not feasi- 
ble to calculate Dn for values of n as large as may be possible with 
the method of reference (l). 
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APPENDIX C 

ALTEHlfATIVE DRAG F(»MULA 


In reference It It is sliown that hodles of revolution exist ’which 
involve only a small numher of values of n in the Fourier series expan- 
sion of the gradient of the hody area dlstrihutlon. Similarly, planar 
configurations exist which Involve only a small, numher of values of n 
in equation (3) when t in equations (l) and (B3) is allowed to vary 
with polar angle cp and the speed parameter p. One such, a wing of 
elliptic plan form and circular-arc secticms, was discussed in refer- 
ence 6 and was shown to have min i mum drag for given volume with elliptic 
plan form. 

In this Appendix the drags of a series of wings of elliptic plan form 
are derived and the drag of an arbitrary planar configuration is expressed 
in teicms of the elliptic-wing drags for the purpose of including the pre- 
dominant effects of spanwise extension of the wing in a small number of 
terms of a series. 

The double moments (Mpk) of an arbitrary planar configuration can be 
expressed in terms of double moments associated with an elliptic plan form 
which encloses the arbitrary plan form. Then the drag of the arbitrary 
configuration will be equal to the drag of the corresponding combination 
of elliptic -wing moments. 

Equations (B 3 ) and (B 5 ) can be combined to obtain 
n-2 n-p-2 

^ " Z Z (f) 

p=o k=o 


Taking I to be the length of the equivalent body of revolution for an 
elliptic plan form of span bj. and maximum chord ci yields 


= ci^ + (pbi COB cp)^ 


or 



1 + 


COB 




COS 


:aHFlDENTIAL 
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With the definition 


Cl = 
I 


cos a 


(C2) 


it is found that 





slnPot C08^+% 


so that 


An 


n-2 


I 

p=o 


I 

k=n 


ehpk M^ic 


slnPa cos^"^^ 


(C 3 ) 


where 



The dimensionless double moments^ (M’plj) of the elliptic Ipu p of 
reference 6 can he found hy setting ^ = 0 for n ^ 2 and A^ = cos^a 
in equation (C3). Similarly, other wings of elliptic plan form can he 
defined hy setting Ajx = 0 for n ni and = cos^a sinPia. Label this 
series of_^wlngs with the nvmhers ni and pi. The corresponding double 
moments Mn^p^^pk of each such wing can he found for even values of n^ 
and pi from the relation 


n -2 


n-p-2 


I I 


gnpk MuiPipksinPa cos^a = d^Q^sinPia 


p=o 

p even 


k=o 

even 


values of 


for all 
a and n 


(C 5 ) 


where 


^nnx 


1 for n = ni 
zero otherwise 


(C6) 


Equation (C5) can only he satisfied for values of Pi ^ ni - 2 and It eaTi 
he seen that the Mnipipk's are zero for p < pi or k < ni - p - 2. In 
the following discussion all integers have even values only. 
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To obtain a drag fornrula, tlie ^M'pk's of the configuration are 
expressed as a combination of the Mn^^p^pk'B of the form 

p p+k+2 

M*pk = E I ^HiPx ^niPipk C^T) 

Pi=o ni=pi+a 

Substitution of this in equation (C3) yields 
n -2 n-p -2 p p+fc+a 

A„.=os% IEEE Snpk Knipi Muj^p^pkain^a cos^a 
p=o k=o Pi=o niPi +2 

or interchange of the order of svmnnations yields 


n -2 


n -2 n-p -2 


An = cos^o, V y ^XPl V y gnpk MniPi^pksinPa cos^ 

J—s t i-J J—i 


Pi=o ni=pi+2 p=Pi k=ni-p-2 

By the use of equation (C5) this becomes 


n -2 n 


An = cos^a ^^xPl j_Bin^^ct 


Px=o ni=pi+2 


Upon completion of the summation with respect to n^ the expression 

n-2 

An = cos^o, ^ KnpiSin^^o, (c8) 

Px=o 

is obtained and can be used in the place of equation (C3) iu the drag 
formula if the found. The fi^p^’s will be linear func- 

tions of the configuration Mpj^'s of the form 


Px n-P2-2 


Ki 


■npr 


= ^ ^ ^npipaka 


(C9) 


P 3 =o k„=o 
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The quantities Kup^ can he found from the terms Mpjj. of the con- 
figuration If the ^niPiPoko's can he found. From equations (Cl) and 
(c8) it is found that ' 


□ -Z n-p- 

I I 


p=o k=o 


n-a 


gnpk WPpkslnPa cos^a = ^ Knp 3 ^sinl*ia 


Pi=o 


for all a and all fimctions M’pt. Substituting equation (C9) In the 
right side of the above equation yields 


n -2 n-p -2 


Pi n^-2 


p=o ki=o 


X X Mpksin^'a cos^ ^ X X X M’pjj,sinP^a 

Pi=o p=o k=pi-p 
n-g n-p -2 k+y 

i5’ 

p=0 kiro Pl=P 


* * ~ ** • 

" X X ^ fnpipksinPia 


or 


n -2 n-p -2 


p=o k=o 


k+p 

X X Tsnpk sln% cos^ - ^ fnpipksln^ia^ = 0 

- > pi=p 


for all a and all functions Mpk. Thus 

k+p 


Snpk®^^^® COSTCO — ^ ‘ Tnpj^pjf sin^^i 


a for all a 


Pi=P 


or expansion of cos^ by the binomial theorem yields 
p+k , k+p 

Pi-P ^ 

^pk 2^ ( “l) ^ ^ Pi-P sln^^ = ^ ^npipk®^^^^ 

Pi=P s P:^ 

for all a so that 


i^cTOTnewriAiiT:: 

* _ J ■ c__. ^ 
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Pi-P t 

^npipk = (-1) ^ Pi-p ^pk 

2 

Substituting equations (l2) and (b 4) in this equation yields 


(CIO) 


fnpipk = 


n-pi-k-2 

(-1) = 


^m+p+k^ . 


2P^k (k/2) i 

pikl 


(Cll) 


for even values of tlie Integers. 

For odd values of Ui and even pi the MnxPiPls's are defined "by 

n-s n-p -2 

I I 8npk ^ipipfcsiii^a- cos^"^ = dnniBlnPia (C12) 

p=o k=i. 
p even 


ari^ hy a similar process 


n-3 


Aji = cos^a KnpsinPa. 


p=o 


Pi n-Pa-2 


Ki 


■npi = 


X X ^^PiPaka ^'pska 


P2=0 k2=Pi-p2+l 
P 2 even kg odd 


(C13) 


(CU) 


and 


Pi-P lEZ± 

fnpipk = (-1) ^ Px-p gnpk 


(C15) 


for odd values of n and k. 


.imoiaiiiiaitMKfJ 
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Since tbe dimensionless drag components are given "by 




integrals of the tsrpe 


1- r 

2n J 


art 


cos qp^ 




sl2+* 
s 


dfp. 


1. r^" 

2 « Jr. 


COS 

Vci 


J 


1 + cos 

\°1 


"T 


lit® 

a 


d£p 


are encountered and have been evaluated hy means of a method of residue 
Integration and differentiation vith respect to parameters. 

In summary, the drag of an arhltrary closed planar configuration with 
spanwlse symmetry can he evaluated hy means of the following formulas: 


D = itq ( t 


CO 

in 

n«=a 


(Cl6) 


Dn = 


an-* 




J, I. 




J=o 

j even 


an -6 


J=o 

j even 


for even n 




(C17) 


= 


Jl. 


'm 


for even j 


(Cl8) 




(EHTIAL . 
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(C19) 





(^J 



3 

(j+2)(j+4) V'^V. 


'•(e) 

2" 

j+5 

2 


(C20) 



p=o 

p even 


(C21) 


I 

P 2=0 

Pg even 


Knp 


P 

z 

P 2=0 

p even 


n-ps-2 


^ ^nppak 


even n 


k2=p-Ps 
ka even 


n-p 2-2 


^ fnppafc; 


■2 ^Pakg 


for- odd n 


ka=p-Pa+^ 
kg odd 


(C22) 






In equation (C24) t(x,y) Is the configuration thickness dlstrlhutlon, 
hi and ci are the span and maxlm-um chord, respectively, of an ellipse 
which completely encloses the configuration, h is the wing span, and c 
is the streamwlse length of the wing. 
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The first five values of Dn . are given "by 
D 2 = 2(P?oo) ^ lo 

fe) 


5s = 3(lfM‘oi)^ Qo 


D 4 , = 4 


(12M'o2 - ^ 00 )^ Ic 


fe) 


+ (12M'o 2 - ^’oo)(l^’ 


20 


l^’oa) I 2 + I (12^20 - 1 ^ 02 )^ I 4 fe)] 

5s = ?[( 3^*03 - ISM'or)^ Qo(^^^+ (S^’os ~ ISM’oi.) (90^2^ - 

32M'o 3) Qa (^) + § (9^2X - 3^‘o 3>^ Q* fe)] 

|(80M»o^ - 48 M*o 2 + 3M*oo)^ To + (SOM’oa " ^’02 ■» 


De = 6 


3M'oo) (480M'g2 - i6cm'o^ + 1(^02 - l<^* 2 o) I 2 


XC25) 


1 [2(80M'o4 •' ^'02 + 3 M’oo) (8CM04 - 1+80M'22 + 80M«4 o) + 

8 

( 480 M *22 “ 160 M'o 4 + Jt8M'Q2 — ^iSM* 2 o )^3 ^4 

^ (i{ac»l'22 - 160M^4 + 4aM'o2 “ WM*2 o) (8CM’o4 “ kSC^zz + 
Sdin^o) le ^ (8 QS'o 4 - i*80M»22 + SCM'zo)^ Is (^)} 


The convergence with respect to n is "best when the smallest possi- 
ble ellipse is used. The theoretical drag cxirves in sketch (e) were cal- 
culated hy dividing the configuration into parts of short, long, and 
Intermediate length so that smaller ellipses could he used for the shorter 
parts. The drags of these three parts and their interactions were calcu- 
lated using values of n up to 6 in the formulas of this Appendix. The 
body moments in the region of the body cutouts were found to be important 
and were taken into account using a quasl-cyllndrlcal approximation to 
find an equivalent planar system neglecting Induced camber effects. 
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